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This is a brief, annual status report which includes the theoretical achievement and the 
experimental setup during the report period of April 1, 1991 to March 31, 1992. 

The theoretical achievement is based on our paper entitled "Rewetting Theory and the 
Dryout Heat Flux of Smooth and Grooved Plates with a Uniform Heating" (see Appendix A) 
which has been accepted in the "Heat Pipes and Thermosyphons" session at the 1992 ASME 
Winter Annual Meeting. Although several investigations have been made to determine the 
rewetting characteristics of liquid films on heated rods, tubes and flat plates, no solutions are 
yet available to describe the rewetting process of a hot plate subjected to a uniform heating. In 
our paper, a model is presented to analyze the rewetting process of such plates with and without 
grooves. Approximate analytical solutions are presented for the prediction of the rewetting 
velocity and the transient temperature profiles of the plates. It is shown that the present 
rewetting velocity solution reduces correctly to the existing solution for the rewetting of an 
initially hot isothermal plate without heating from beneath the plate. Numerical solutions have 
also been obtained to validate the analytical solutions. In Appendix A, the successful prediction 
of the rewetting curve by the approximate analytical solution for the grooved (or smooth) plate 
with a uniform heating is illustrated. It shows that the approximate closed form solution is in 
reasonably good agreement with the numerical solution. Furthermore, a simple method is 
presented in Appendix A to predict the dryout heat flux of a liquid film flowing over a heated 
smooth or grooved plate. The results of the prediction are found to be in reasonable agreement 
with the existing experimental data. 

During this first annual period, an experimental system has been set up, design of which 
has been reported in our semiannual status report on October 18, 1991, which includes (1) 
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microfin surface geometry of the grooved plate; (2) the selection of the test fluid and plate 
material; (3) surface temperature measurement mechanism; (4) liquid flow rate control and 
measurement design; (5) heat source of constant temperature/variable heat flux system. The 
complete experimental system including the data acquisition system was being built and 
assembled during the second half of the year. Some shake-up experiments to measure rewetting 
velocity and the surface temperature distribution of the plate will be conducted. Figure 1 shows 
the overall experimental system while Figure 2 shows the test section. Some experimental 
results will be reported in the next semiannual status report. 



Fig. 1. Heat Pipe Rewetting Experimental System 
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Fig. 2. Test Section 
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APPENDIX 


A reprint of the paper entitled, "Rewetting Theory and the Dryout Heat Flux of Smooth 
and Grooved Plates With a Uniform Heating," to be presented at 1992 ASME Winter Annual 
Meeting. 


To be presented at 1992 ASME Winter Annual Meeting, Anaheim, CA. 


Nov. 8-13, 1992. 
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REWETTING THEORY AND THE DRYOUT HEAT FLUX OF SMOOTH 

And grooved plates with a uniform heating 

S. H. Chan and W. Zhang 
Department of Mechanical Engineering 
University of Wisconsin— Milwaukee 
P. O; Box 784 
Milwaukee, Wisconsin 53201 


ABSTRACT 

The evaporation and condensation of thin liquid films are 
of significant importance in a wide variety of problems ranging 
from specific applications in the heat pipe field to more general 
ones in chemical, nuclear and petrochemical industries. Although 
several investigations have been conducted to determine the 
rewetting characteristics of liquid films on heated rods, tubes and 
flat plates, no solutions are yet available to describe the rewetting 
process of a hot plate subjected to a uniform heating. A model is 
presented to analyze the rewetting process of such plates with and 
without grooves. Approximate analytical solutions are presented 
for the prediction of the rewetting velocity and the transient 
temperature profiles of the plates. It is shown that the present 
rewetting velocity solution reduces correctly to the existing 
solution for the rewetting of an initially hot isothermal plate 
without heating from beneath the plate. Numerical solutions have 
also been obtained to validate the analytical solutions. Finally, a 
simple method is presented to predict the dryout heat flux of a 
liquid film flowing over a heated smooth or grooved plate. The 
results of the prediction are found to be in reasonable agreement 
with the existing experimental data. 


NOMENCLATURE 

A dimensionless heat source 

coefficient defined in eq. (65) 

B Biot number with respect to the convective heat transfer 

coefficient 

C p thermal capacitance (J/Kg-°C) 

h surface convective heat transfer coefficient (W/m z -°C) 

h average boiling heat transfer coefficient (W/m 2 -°C) 

h* convective coefficient between plate and environmental 

gas (W/m 2 -°C) 


K thermal conductivity (W/m-°C) 

L length of the plate (m) 

Lj length of wet region (m) 

Lj length of dry region (m) 

N grooved geometric coefficient 

n exponent defined in eq. (65) 

P dimensionless rewetting velocity 

q uniform heat flux (W/m 2 ) 

^Ichf critical heat flux (CHF) (WYm 2 ) 

incipient boiling heat flux (W/m 2 ) 
q^ maximum heat flux (W/m 2 ) 

s plate thickness for smooth (sj) or grooved (s { - i ,) 

plates (m) 

Sj plate thickness (m) 

s 2 liquid film thickness (m) 

T temperature (°C) 

T a environmental temperature (°C) 

T f inlet liquid temperature (°C) 

Tj incipient boiling temperature (°C) 

T 0 Leidenfrost temperature (°C) 

T s saturation temperature (°C) 

Tj initial hot surface temperature (°C) 

t time (sec) 

U r rewetting front velocity (m/s) 

x length of liquid film in moving front coordinate (m) 

x' length of liquid film in stationary system (m) 


Greek symbols 

constant defined in eq. (23) 

<* 2 constant defined in eq. (45) 

constant defined in eq. (23) 

0 2 constant defined in eq. (45) 

rj dimensionless length coordinate with respect to x 

i?' dimensionless length coordinate with respect to x' 

9 dimensionless temperature 

0 h dimensionless temperature of the transient part 

9 S dimensionless temperature of the steady-state part 

r dimensionless time 


initially wetted by a liquid film. A means is provided for the 
prediction of the maximum heat input which results in the dryout 
of the liquid on the plate. The result of the prediction is compared 
with the experimental data of Grimley (1988). 

REWETTING MODEL AND SOLUTIONS 

Consider the rewetting process of a hot plate initially at a 
uniform temperature Tj with no liquid on the plate as shown in 
Fig, 1, the plate is heated below by a uniform heat flux and is 
quenched by a liquid advancing along the direction of the grooves 
on the top surface of the plate. In order to simplify the complexity 
of the physical phenomena of the rewetting process, we consider 
first a heated smooth surface plate with no parallel grooves. 


P density (Kg/m 3 ) 

Subscripts 

d dry region 

g grooved plate 

L total length of the plate 

L| length of wet region 

L 2 length of dry region 

s smooth plate 

w wet region 


INTRODUCTION 

The rewetting process is a conjugated heat transfer problem 
involving interactions between a solid wall and flowing fluids. 
The process for rewetting of a grooved plate with a uniform 
heating is complicated, as the rewetting velocity varies with time, 
physical geometry of the grooves, plate properties, fluid properties 
and the applied heat flux. . Although several investigations 
(Yamanouchi, 1968; Thompson, 1972; Duffey, 1973; Sun, 1974; 
Alario, 1983; Grimley, 1988; Stroes, 1990; Femg, 1991; Peng 
and Peterson, 1991) have been made to determine the rewetting 
characteristics of liquid films on heated rods, tubes and flat plates, 
none has yet presented the solution for the rewetting process of a 
heated plate with a smooth or a grooved surface subjected to a 
uniform heat flux. It is noted that the surface with small grooves 
has received increasing attention as it has many practical 
applications. For instance, the microgrooved surface is employed 
most often to enhance heat transfer (Grimley, 1988). 
Microgrooves are also useful for replacing the wicking material in 
heat pipes. In fact, it is used in the innovative monogroove heat 
pipe design for the thermal radiators of the space station (Alario 
et al., 1983). The recent attempt to solve for the rewetting 
process of a plate subjected to a uniform heating by Peng and 
Peterson (1991) has apparently encountered a difficulty as their 
transient temperature profile of the heated plate was incompatible 
with their boundary conditions. One of the objectives here is to 
present a physical model suitable for the rewetting analysis of a 
plate heated by a constant heat flux from below. Solutions for the 
transient temperature profile and the rewetting velocity are 
presented. The incompatibility problem is also resolved. Another 
objective is to analyze the diyout limit of a grooved surface 


It Smooth Surface Plate 

The rewetting process of a hot dry plate is sketched in Fig. 
1. A liquid film from a liquid reservoir, driven by its surface 
tension, is to advance along the hot plate. Similar to prior studies 
by Yamanouchi (1968), Duffey and Porthouse (1973), and Sun et 
al. (1974), the initial temperature of the plate is assumed to be 
higher than Leidenfrost temperature T 0 such that the rewetting 
process is assumed to be conduction-controlled. Also, an averaged 
heat transfer coefficient is assumed in the wet region to remove the 
heat from the thin plate to the liquid film, no heat loss to the 
environment is assumed in the dry region, the plate at the rewet 
front is assumed to remain at a constant Leidenfrost temperature, 
and the plate is thin enough that the one-dimensional rewetting 
model can be invoked. It is therefore proposed to solve the 
following governing equation (see Appendix A) on a Lagrangian 
coordinate moving with the rewetting front, 


99 m a 2 fl + p de 
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where, B is the Biot number; A, P, 0 ( 17 , 7 ), 77 and 7 are the 
dimensionless heat source, rewetting velocity, temperature, length 
and time respectively. In the wet region (-Lj ^ x ^ 0), 
h=constant?*0 while in the dry region (0 < x < I^), h=B~0. 

The above governing equation is different from the equation 
in ail prior rewetting models (Yamanouchi, 1968; Duffey, 1973; 
Sun, 1974; Peng, 1991 et al.) in that an extra term, namely, the 
transient term on the left-hand side of the equation, is added. The 
addition of this transient term is essential in analyzing the 
rewetting process of the heated plate subjected to a heat flux 
conduction as shown in Fig. I. This is because the temperature 
profile of the plate in the Lagrangian coordinate (x, y) is no longer 
invariant with time as in the case of rewetting an initially hot 
isothermal plate without heating from beneath the plate. Without 


it, the incompatibility difficulty as noted above will arise and the 
final solution can not be expected to satisfy the transient boundary 
condition at x = <». The addition of the transient term, however, 
renders some mathematical complications in the solution of the 
re wetting velocity. 

The initial condition is 


e (17,0) 


- T s 
T 0 - T, 


while the boundary conditions are, 

e (~v Ll ,r) = 0 
e (o,r) = l 


8 0?V r ) = T~ - r ‘ e ' + At 


where i 7 L1 and ij L2 are defined as 

_ L l - L 2 
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and T l is the dry plate temperature at ij L2 . The condition given in 
eq. ( 6 ) implies that the plate is long enough such that, within the 
time period of interest, the plate temperature at t ? L2 is not 
thermally affected by the rewetting process. Consequently, it can 
be expressed as 


Tl = Ti + 


¥ 


p C p s i 


( 8 ) 


which is used in eq. ( 6 ). 

In the wet region (-77 L1 ^ rj £ 0), the mathematical model 
of this problem is given as 


dr q v 2 drj 
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8 (V,r) = 6 S ' W (,) + d Kw ( v ,r) (13) 

where 0 S W is the solution for the steady-state part of the problem, 

(14) 
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8 s,w (~Vl,) = 0 (16) 

while 0 h w is the solution to the transient part of the problem 
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The solution of the steady-state problem is given as 
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To solve for the exact, analytical solution of the above 
appears to be difficult, if not impossible. Therefore an 
approximate, analytical solution is sought. This is made possible 
by treating the Peclet number, P, as a constant value in the 
mathematical deliberation in order to achieve a closed form 
solution. This approximation appears to be reasonable as the 
rewetting velocity tends to reach a quasi-steady state after an initial 
period when the liquid film is brought in contact with the hot 
plate. The numerically exact solution will also be presented later 
to check the accuracy of the closed form solution. Accordingly, 
the solution is split into two parts, 


r _ -p - / p 2 7 
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To solve for the transient solution, the following 
transformation is introduced, 


8 h.w (.V,r) = e a ' T ' * B,t p(tj , 7 ) 


( 22 ) 



where 


0 (V,r) = e sw (v) + e h w (r/,r) 


«i = ~P /2 , fit = -B - P 2 1 4 (23) 


Then eqs. (17) to (20) become 
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The solution of the above can be readily obtained from 
Carslaw (1959) as follows: 


* -nVr/r/2 

' (V,r) = £ a sin (-Oil) e L ‘ (28) 
«=i 1 % 


where 

[° /,OOsin(-dlZ)^ (29) 

J '% V, 


where a,,! and fj(V) are given by eqs. (29) and (30) respectively, 
a x and are given by eq. (23). 

Before proceeding to the solution for the dry region, it is 
noted that the above steady-state solution is identical to the solution 
presented by previous workers (Sun et al., 1974) when A = 0 (or 

q “ 0). 

Similarly, the problem in the dry region (0 <^rj <_ rj L2 ) is 
described by 


do a 2 e + p de 

dr dr] 


+ A 


(32) 


6 (,7,0) = e x 


(33) 


/,OO=e'“' V {0 1 - 


(l-d) e ^ + d 

B_ B 

B _ ‘Vta, 


1-d- 


r,V 


(30) 


(1 + d )e + d 
B J B 


e ~ r *h - e ~ r *h 


Therefore, the combined solution is 
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which is likewise split up into the steady-state and transient 
solution, 
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The solution of the steady-state problem is 
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A similar transformation is introduced 

0 = e aj ’ ! ^ T «( r), t) 
where: a 2 = -pi 2, 0 2 = -p z /4 

to reduce eqs. (39) to (42) to the form 
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which can be solved by the method of linear superposition. Hie, 
the transient temperature profile in the dry region can be readily 
found as 


oo 

' 

E Sin 

nmr) 

n = l 

V 


JL (At) e ♦ ± £ (- 1)» 5in<? " ,? 


’It, «-i 


,-A^mLi . 6 l + p^-l 


e d4> 2 (k)\ * I - 


* A 

i + -p v h~ 
e'^-l 


e \ + j % - 1 

e'^-1 


- -17 

n ' 


where /?„ = nx/i| L2 and 


«„2 « — / 2 (w) sin ^ dw 

% i % 

At the rewetting front, the conductive heat flux is 
continuous, i.e., 

— \ - (5D 

dr > L- dr > Uo* 

Upon the substitution of eqs. (31) and (50) into eq. (51), the 
dimensionless rewetting velocity P can be determined by the 
expression 
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It is of interest to examine the limiting solution of the 
above for the case that has been investigated by Yamanouchi 
(1968), namely, the rewetting of an infinitely long, hot and 
isothermal plate without any heating. In this case, A=0 (i.e., 
q=0) and then by setting 7j L1 = 17^2 00 , t -* 00 # the above 

solution reduces to 
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Combining Eqs. (54), (55) and (56) yields 


( 56 ) 


which is exactly the same as the well known Yamanouchi’ s 
solution (1968). 

Due to the absence of data on rewetting velocity on a hot 
plate heated by a uniform heat flux beneath the plate, the solution 
given by eq. (52) can not be compared with experimental data. 
However, experimental data are available on a hot plate without 
heating from below. Yamanouchi (1968), for example, has 
confirmed reasonable agreement between the limiting solution 
given by eq. (53) and his data. The general solution presented 
above for a smooth plate will be used in the following section to 
yield the solution for the plate with axial grooves. 

2 . Grooved Plate 

The rewetting model of the grooved plate is based on that 
of the smooth surface plate. The coolant is driven by the wicking 
(surface tension) effect of microfms and is assumed, without loss 
of generality, to fill up the grooves as shown in Fig. 2. 

At the level of y = 0, 



(2f 2 +f 3 ) = h(T - T s )(2( x +f 3 ) 

(54) 
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where the grooved geometric coefficient is defined as 


N = (2t x + £ 3 )/(2 £ 2 + £ 3 ) 

The result given by eq. (57) suggests two useful 
simplifications. First, the factor (Nh) is the equivalent convective 
heat transfer coefficient of the grooved plate and can be 
approximated by that of a thin liquid film on a smooth plate 
multiplied by a factor of N, namely, 

[ heat transfer coefficient 
of a grooved plate 

where N is more generally defined as 

^ _ fhe wetted perimeter 
width of the cross section 


heat transfer coefficient 
of a smooth plate 


+ h a (T-TJ- 21 2 

where, h is the convective coefficient of a smooth surface, T a is 
the environmental, temperature, h a is the convective coefficient 
between the plate surface and the environmental gas above ( 1 ^ = 
0). At the level of y =» s t - £ lf 


— \ = 1 

dy [y - k 
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Since the plate is thin, 


Second, the governing equation for the grooved plate with a 
uniform heating remains unchanged provided h is replaced by (Nh) 
and the dimensionless variables are properly scaled as follows: 

P = U r p C p (s, - i x )!K ; B - Nh (s x - i x )f K 

A -q (s Y ~ t x )l[K(T a - T s )); v =x/(s { - £,) 

^ =L,/(Ji - f,); % “V(*l - <,) 

T = t/tet - efpCp/K] (59) 

where h is the convective coefficient of the smooth plate with a 
uniform heating. Therefore, the solution of the grooved plate with 
a uniform heating in the wet region is the same as eq. (31) with 
the above modification and the solution for the dry region is 
identical to eq. (50). 

Numerical solutions have also been obtained by solving for the 
original governing equation (see Appendix A, eq. (A7)) fixed to 
the nonmoving coordinates (x', t), 
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discussed above. However, they did report interesting data on the 
maximum heat flux (the critical heat flux, CHF) that the heater 
could supply to the plate without causing the diyout of the flowing 
liquid film. Since the dryout and the subsequent rewet of a heated 
surface are an integrated problem in heat pipe applications, it is 
desirable to be able to explain or predict the heat flux condition 
that leads to dryout of the plate. To achieve this objective, a 
simple method is presented next and comparisons will be made 
with the reported data. 


PREDICTION OF THE MAXIMUM DRYOUT HEAT FLUX 


where 


V ■ x' Is, i) r L s l Is 

and s = Sj and (s x - t x ) for the smooth and grooved plates 
respectively. The rewetting front location, is determined 
from 0(V L1 , r) = 1, The differential equation is discretized by a 
standard finite-difference approach. Figure 3 illustrates the 
computed wall temperature of the grooved (or smooth) plate versus 
the plate length, the so-called rewetting temperature curve, at 
various times when the liquid film is FC-72, which is Case 1 of 
Fig. 4. These numerical values are selected from the experimental 
condition of Fig. 5 of Grimley et al. (1988), namely, for FC-72 
fluid on the smooth surface, T { - T s = 11°C, T 0 - T 8 = 17.5°C, 
T s = 56°C, Tj = 100°C from which 6 t is calculated. Other 
dimensionless parameters A, B and q L are estimated from the 
following values. The needed heat transfer coefficient h =2779.49 
W/m 2o C is calculated from their boiling curve by eq. (65). For the 
grooved copperplate, C = 383.1 J/kg°C, p = 8954 kg/m 3 , K - 
386 W/m°C. When the Liquid fills up the grooves as shown in Fig. 
2, the grooved geometric coefficient is N = 1.75, which is based 
on their geometric dimension of t x = 0.5mm, t 2 = 0.2mm, f 3 = 
0.4mm, Sj = 6mm, and q = 24000 and 100000 W/m 2 for case 1 
and 2 respectively. They also reported dryout heat flux data which 
will be used in the later comparison. From Fig. 3, it is clearly 
shown that the rewetting temperature profiles are transient in 
nature, even in the coordinate frame moving with the rewetting 
front. 

The successful prediction of the rewetting curve by the 
approximate analytical solution, eq. (52), for the grooved (or 
smooth) plate with a uniform heating is illustrated in Fig. 4. It 
shows that the approximate closed form solution is in reasonably 
good agreement with the numerical solution. This may appear to 
be somewhat of a surprise in view of the seemingly inconsistent 
approximation of treating P (or Ur) as constant. As mentioned 
above, such an approximation was necessarily made in the 
mathematical manipulation to achieve an approximate closed form 
solution. It is equivalent to neglecting higher order terms attributed 
to the transient components of P. Figure 4 shows that such an 
approximation does yield good results as is expected because all 
rewetting velocities tend to level off quickly. In fact, beyond 
t= 5.208 (or t=1.4 sec), the difference between the closed form 
and the numerical solutions is also indiscernible. 

In a recent study, Grimley et al. (1988) conducted 
experiments to investigate the enhancement of convective boiling 
heat transfer by grooves on a plate heated from beneath. 
Unfortunately no rewetting velocity data were reported which 
could otherwise be useful to check the validity of the solutions 


Under consideration is a smooth or a grooved plate initially 
covered by a thin flowing liquid film at a temperature T s . The 
plate is subjected to a uniform heating. It is of interest to predict 
the maximum heat flux that triggers the dryout of the film. 

A simple method based on the above rewetting concept is 
now extended to provide a means of estimating this maximum 
dryout heat flux. For a flowing liquid film over a plate heated by 
a heat flux which exceeds the maximum heat removal capability by 
convection and boiling, the liquid will cease to advance and begin 
to recede. Thus dryout will occur. On this physical premise, the 
maximum rate of the heat removal is given by 

= * (T 0 ,s - D 

where, following the work of Howard et al. (1975), the average 
convective boiling heat transfer coefficient of the liquid film is 
estimated from 


h 




(65) 


in which T t is the plate temperature at the onset of boiling, and 
T 0 8 is the smooth plate Leidenfrost temperature of the rewetting 
front and is approximated by the plate temperature at the CHF 
location. Q b is the boiling curve of the liquid film over the plate. 
As an approximation, a form Q b =aj(T - Tg)", which fits the 
boiling curve, can be used. In the event that the boiling curve of 
the flowing liquid film is unavailable, the pool boiling curve could 
be used as the first order approximation (Howard et al. 1975). 

In the case of the grooved plate, the same analogy developed 
above is adopted here. The maximum dryout heat flux is estimated 
from 


a = N 


K * (To. - T s ) 


( 66 ) 


where T 0 g is the Leidenfrost temperature of the grooved plate, 
also approximated by its CHF temperature. Based on the 
experimental data of Grimley et al. (1988) for a fluorocarbon (FC- 
72) liquid film falling over heated smooth and grooved plates, T 0 g 
was found to be slightly lower than T 0 s . Thus, if T 0 g of the 
grooved plate is unavailable due to the lack of the boiling curve 
for the grooved plate, one may attempt to approximate T 0 from 
the smooth surface data T 0 s . Then the predicted maximum heat 
flux may be slightly overestimated, namely, 

<7ma* .g < N • Tl ■ ( T 0,s ~ D 


As an application to show the feasibility of the above 
method, the experimental conditions and geometries of Grimley et 



al. (1988) are used. In their experiments, the liquid film 
completely covered the smooth and grooved plates, ^=0,5 mm, 
*2=0.2 mm and * 3 =0.4 mm, such that N=2.25 (from eq. (58)). 
(Ti - Tg) and (T 0 - Tg) are taken from their boiling curves. The 
average boiling heat transfer coefficient is calculated from their 
smooth surface boiling curve and eq. (65). The two correlation 
constants aj and n are determined by arbitrarily collocating the 
boiling curve of the smooth surface plate at two locations, T { and 
T 0 , where the boiling heat fluxes are designated by and q CHF , 
respectively. Table 1 shows the comparison between the predicted 
maximum dryout heat flux, q max , using eqs. (64) and (66), and the 
reported dryout data q CHF for both types of plates. The agreement 
is satisfactory particularly in view of the simplicity of the method 
proposed for the grooved plate. The same agreement is shown in 
Table 2 when the mean inlet velocity of the falling film is 
increased from 0.5 m/s to 1.0 m/s. 


CONCLUSIONS 

The rewetting process of a smooth surface plate subjected 
to a uniform heating has been investigated. A proper governing 
transient heat conduction with a convective boiling condition has 
been solved to yield an approximate closed form solution for the 
plate temperature profiles in the wet and dry regions of the plate. 
From the temperature profiles, an approximate closed form 
solution for the rewetting velocity over the heated plate has been 
obtained. Numerical solutions have also been presented to check 
the validity of the closed form solution. The closed form rewetting 
velocity was found to be in good agreement with that of the 
numerical solution. It is shown that in a limiting condition the 
present rewetting velocity solution reduces correctly to the existing 
solution for the rewetting of a hot, isothermal plate without 
heating. However, contrary to the case without heating, the 
rewetting process on the plate with a uniform heating is found to 
be transient (time variant) even on the coordinate frame moving 
with the rewetting front. The rewetting velocity is found to be 
much faster initially and then levels off later. 

A method to address the rewetting process of the grooved 
plate based on the smooth plate rewetting model has been 
developed. It has been shown that, by properly defined scalings, 
the solution for the smooth plate can be made to be applicable for 
the grooved plate. 

Finally, the dryout of a liquid film over a heated plate has 
been investigated. A simple method has been proposed to predict 
the dryout critical heat flux of the smooth and grooved plates. The 
results of the prediction were compared and found in reasonable 
agreement with the existing experimental data. 
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APPENDIX A: CO-ORDINATE TRANSFORMATION 

The heat conduction equation within a smooth plate as 
shown in Fig. 1 is 



d 2 T + d^T_ 
3x' 2 dy' 2 


(Al) 


For a coordinate system x-o-y moving with the rewetting front, 


! u r (t)dt 


(A2) 
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eq. (Al) can be written as 


PC, 


’fi 

I] -U r 

' dT 


= K 

9 2 T + 9 2 T 

‘ 

)i U 

dx 

t 


dx 2 9 fy 


(A5) 


For a thin plate subjected to heating from below, it further reduces 
to 



Fig. 2. Microfin Structure and Cross-section Shape of Grooved 
Plate 


far] „ _k_ £r + u dT 

[to lx PC„ dx 2 r dx 

or in dimensionless form 


h 

P c ps i 


(T~T t ) 


♦ 1 
PC/i 


»-*+pZ-B9*A (A6) 

dr dr} 2 dr, 

where the dimensionless variables and parameters are given by eq 

( 2 ). 

For the purposes of comparison and numerical computation, 
the above is written in untransformed coordinates (x', t) as 


98 

dr 


d 2 0 

3v 2 


- B8 + A 


(A7) 


where 

v' m x'/s, (A8) 



Fig. 3. Numerical Solution for Grooved Plate 
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Fig. 1. Schematic of the Rewetting Film on the Heated Plate 



Fig. 4. Rewetting Velocity Comparison, <91 —2.5143, 
ij, = 14.5455, B=0. 0693 (wet region), A=0. 01954 (Case 
1), A=0. 08142 (Case 2). 






Table 1. Comparison Between Predicted and Measured (Grimley, 1988, Fig. 8,9) 
Maximum Dryout Heat Flux 



Table 2. Comparison Between Predicted and Measured (Grimley, 1988, Fig. 5) 
Maximum Dryout Heat Flux 
















































